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Abstract 

We describe a calculation of the 0{a 2 s ) QCD corrections to the fully differential cross section for W and 
Z boson production in hadronic collisions. The result is fully realistic in that it includes spin correlations, 
finite width effects, 7 — Z interference and allows for the application of arbitrary cuts on the leptonic 
decay products of the W and Z. We have implemented this calculation into a numerical program. We 
demonstrate the use of this code by presenting phenomenological results for several future LHC analyses 
and recent Tevatron measurements, including the W cross section in the forward rapidity region and the 
central over forward cross section ratio. 



1 



I. INTRODUCTION 



The study of the electroweak gauge bosons W and Z is an important part of the physics pro- 
grams at the Tevatron and the LHC. Their large production rates and clean experimental signatures 
facilitate several important measurements, such as the determination of the electroweak parame- 
ters M w and sin 2 6^, and the extraction of the parton distribution functions of the proton [1]. The 
idea to use W and Z production to monitor hadron collider luminosities has begun to be studied 
at the Tevatron and will continue to be investigated at the LHC [2]. During the initial stages of 
LHC running, W and Z bosons will be used to calibrate lepton energy scales, test the uniformity 
of the electromagnetic calorimeter and study the tracker alignment [3]. Searching for deviations 
from Standard Model predictions in di-lepton events with large invariant mass, missing energy, or 
transverse momentum probes extensions of the Standard Model which contain new gauge bosons 
or other exotic resonances. 

The extensive experimental program described above requires accurate theoretical predictions, 
and many calculations describing electroweak gauge boson production are available. For example, 
the fully differential level 0(a s ) next-to-leading order (NLO) QCD [4] and one-loop electroweak 
corrections |0] to W, Z production have been known for many years. The NLO QCD corrections 
are incorporated into the parton-shower event generator MC@NLO [6], which provides a consis- 
tent description of 0(a s ) hard emission effects together with a leading-logarithmic resummation 
of soft and collinear QCD radiation. The j9 T -spectra of W and Z bosons are computed with next- 
to-leading-logarithmic accuracy and are incorporated into the program RESBOS 021]. The 0(a 2 ) 
QCD corrections in the large transverse momentum region are also known [ 8]. However, for many 
applications these theoretical results are insufficient. With an integrated luminosity per experiment 
at the Tevatron exceeding 1 fb _1 and an expected luminosity at the LHC of about 10 fb _1 /year, 
the statistical error on W, Z production is becoming smaller than one percent. This statistical error 
sets the scale for the desired theoretical precision, so that the physics potential may be fully ex- 
ploited. It is then easy to see that for percent level accuracy, next-to-next-to-leading order (NNLO) 
QCD computations for electroweak gauge boson production are required. The available results for 
0(a 2 s ) corrections to inclusive W, Z production B9D and to W, Z rapidity distributions [ 10] confirm 
that the NNLO QCD effects are at the level of a few percent, and that the remaining theoretical 
error after these corrections are included is at the percent level or lower. 

The experimental identification of the W and Z production processes requires cuts on the 
pseudorapidities and transverse momenta of charged leptons and on the missing energy. The 
NNLO QCD results to the inclusive cross sections and to the rapidity distributions cannot be used 
to calculate the effects of these cuts, since they do not contain the spin correlations between the 
leptons and the initial- state partons arising from the spin-one nature of the electroweak gauge 
bosons 

O. The fully differential 0{a 2 s ) corrections with spin correlations included are required 
to model these effects. To illustrate this point further, we discuss three examples. 

1. CDF and DO have recently presented Run II measurements of the W/Z cross section ra- 
tio lfl2ll . With only 72 pb _1 of integrated luminosity, CDF obtained 

<7w x Bi(W —> liA 

R = — - -f- = 10.92 ± 0.15 stat ± 0.14 sys , l = e,fi. (1) 

cr Z / 7 * x Br(z — > l + r) 

The systematic error was estimated to be only 1.5%. The actual measurement was per- 
formed with the cuts p T > 20 GeV and |^| < 1.8 on the charged lepton, and with missing 
energy fi T > 25 GeV. The results were extrapolated to total cross sections using theoreti- 
cally computed acceptances. Since no NNLO QCD calculation capable of modeling these 
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cuts was available at the time of that analysis, the procedure for obtaining acceptances was 
to reweight the PYTHIA rapidity distributions for the electroweak gauge bosons with the 
NNLO QCD computation OH to account for NNLO QCD effects. While it is unlikely that 
this procedure leads to drastically wrong results, a precision of a few percent cannot be 
guaranteed using this technique. 

2. The Tevatron and the LHC can potentially provide stringent constraints on parton distri- 
bution functions (PDFs). The current PDF extractions permit computations of LHC hard- 
scattering cross sections with Q rs 100 GeV to an accuracy of 5% or better. Measurements 
such as the W, Z-charge asymmetries lfl3ll can reduce these errors. These require precision 
predictions through NNLO in QCD for leptonic pseudorapidity distributions with cuts on 
transverse momenta and missing energy. 

3. With the high luminosity of the LHC, precise measurements of electroweak parameters are 
possible. An interesting example is the measurement of the effective electroweak mixing 
angle sin 2 6 W through the forward-backward asymmetry in Z — > l + l~ . A precision of 
2 x 10~ 4 , competitive with the LEP analysis, can be achieved with 100 fb _1 of integrated 
luminosity, provided the following cuts on the leptons can be imposed: |^ e +, e -| < 2.5, 

| V(e + , e") l> 1, p l T > 20 GeV. Given the magnitude of the NLO QCD corrections for this 
set of cuts 11411 . the inclusion of NNLO QCD effects seems mandatory. 

The above discussion illustrates the importance of having a fully differential description of 
electroweak gauge boson production through 0(a 2 s ) in QCD. Unfortunately, such computations 
remain challenging. Their complexity lies in the intricate structure of soft and collinear singular- 
ities that plague individual contributions in QCD perturbation theory. While at NLO a number of 
approaches [15] can be used to isolate and subtract those singularities from complicated matrix 
elements in a process- and observable-independent way, an extension of this approach to NNLO is 
not complete 11611 . We have formulated an alternative technique in a recent series of papers W17I1 . 
The central idea of this method is an automated extraction of infrared singularities from the real 
radiation matrix elements and a numerical cancellation of these divergences with the virtual cor- 
rections. We have previously applied this approach to the computation of the fully differential 
Higgs boson production cross section in gluon fusion and to e + e~ — > 2 jets through NNLO H17I1 . 
We make extensive use of these references in this manuscript. Our goal is to describe a fully real- 
istic calculation of single electroweak gauge boson production at the Tevatron and the LHC. The 
computation is valid through NNLO in perturbative QCD, includes spin correlations, finite widths 
effects, 7 — Z interference and is fully differential. A short version of this paper with initial results 
was presented in 111 811 . 

This manuscript is organized as follows. In the next Section we briefly recall the important 
features of the method and discuss some of the differences between the current calculation and the 
fully differential computation of pp — ■> H — > 77, reported in 11711 . In Section III we demonstrate 
the possible uses of our numerical program by presenting phenomenological results for Tevatron 
and LHC measurements. We conclude in Section IV. 

H. DETAILS OF THE COMPUTATION 

We consider the production of a lepton pair in hadronic collisions, 

h(Pi) + h 2 (P 2 ) -^V + X ^h + l 2 + X, (2) 
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where V = W, Z and i 1)2 represent charged leptons or neutrinos, as appropriate. Within the 
framework of QCD factorization, the cross section for this process is 

da v = ^ / dx 1 dx 2 fj ll {x 1 )fj 2 {x2)dcr ij ^ v+ x{x 1 ,X2), (3) 

where the //* are parton distribution functions that describe the probability to find a parton i with 
momentum xPh in the hadron h. The partonic cross sections der^ are computed perturbatively as 
an expansion in the strong coupling constant a s : 

*°« = <#.v + (*) + og!. v + O(o"). (4) 

The partonic processes i + j — > Z x + Z 2 + -X" that contribute to electroweak gauge boson production 
differ at each order in the perturbative expansion. At leading order in a s , only quark-antiquark 
annihilation channels contribute, while at NLO the (anti)quark-gluon channels also occur. At 
NNLO, the gluon-gluon fusion and (anti)quark-(anti)quark scattering channels also contribute. 
We note that all the relevant channels are catalogued in great detail in Ref. B2D - 

Each partonic process contains several distinct contributions. We describe these components 
using Z-boson production in qq annihilation as an example. At NNLO, this process receives 
three distinct contributions: (i) the two-loop virtual corrections to qq — > Z; (ii) the one-loop 
virtual corrections to qq — ► Z + g; (iii) the tree-level processes qq — ► Z + gg, qq — > Z + qq. 
We refer to these three mechanisms as the double-virtual, the real-virtual, and the double-real 
emission corrections. When computed separately, these three contributions exhibit numerous soft 
and collinear singularities. To produce a physically meaningful result, they should be combined in 
the presence of an infrared-safe measurement function. In addition, the collinear renormalization 
of the parton distribution functions //* is required. 

The double- virtual corrections are the simplest to calculate. They require the two-loop massless 
triangle diagrams that were obtained in [19]. If dimensional regularization is used to regularize 
ultraviolet, soft and collinear singularities, the result is given by a Laurant series in the regular- 
ization parameter e = (d — 4)/2, where d is the dimensionality of space-time. The situation is 
more complex for the real-virtual and double-real corrections because they contain real emission 
matrix elements. These are finite for non-exceptional final-state momenta, but they diverge once 
an emitted parton becomes either soft or collinear to another parton. The challenge in performing 
NNLO computations is to extract the divergences from the real emission matrix elements with- 
out integrating over any kinematic parameter that describes the real emission process. The fully 
differential nature of the computation remains intact only if this can be achieved. 

We have developed a technique to accomplish this in a previous series of papers d We 
describe here the salient features of this method. Consider a double-real emission contribution to 
the production of a Z-boson in qq annihilation, q(pi) + q{jp%) — > Z(p z ) + g{pz) + 9{Pa)- We 
choose a parameterization of the final-state momenta that maps the allowed phase-space onto the 
unit hypercube: 

1 5 

J d d p z d d p 3 d d p 4 5 + (p 2 z - M 2 z )5 + (pl)5 + (pl)5( Pl +p 2 -p 3 -p 4 -p z ) = J JJdA i F({A i }). (5) 

o 1=1 

The function F({Aj}) depends on the details of the parameterization. The invariant masses of 
all particles that participate in the process, such as (p 3 + p 4 ) 2 , (p 3 + p z ) 2 , etc., become functions 
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of the parameters Aj. Soft and collinear singularities in the matrix elements occur when some 
invariant masses reach zero or other exceptional values: (p3 + p^) 2 — > 0, (ps + pz) 2 — > M|, etc. 
Those limits correspond to the edges of phase-space and generally occur when a subset of the Aj 
approaches zero or unity. Two things can happen in these limits. Preferably, the singular limits 
occur in a "factorized" form, and the singularities can be extracted using the simple prescription 
for plus-distributions: 

n=0 

If the singular limits do not factorize but instead appear in an "entangled" form, such as l/(Ai + 
A 2 ), the singularities are disentangled using iterated sector decomposition [20]. We find that at 
NNLO all singular limits can be reduced to one of these two forms. When applied to real emission 
diagrams, this procedure enables us to rewrite the real emission contribution as a Laurant series in 
the regularization parameter e. The coefficients of this series can be integrated numerically over 
phase- space in the presence of arbitrary kinematic constraints. The double- virtual, real-virtual and 
the double-real emission contributions can then be combined and the cancellation of divergences 
for infrared safe observables can be established numerically. A detailed discussion of the method 
with examples of parameterizations used in actual computations can be found in Ref. W17I1 . We 
describe here a few novel aspects of the current calculation. 

• Since electroweak gauge bosons couple to fermions chirally, we must specify our treatment 
of the axial current in d-dimensions. This issue arises from Dirac structures of the form 
TrH[r^75]TrL[r^75], where r^ 1 * 2 ) denote generic products of Dirac matrices and Trn,L 
refer to traces over hadronic and leptonic degrees of freedom, respectively. Unlike in fully 
inclusive computations, these traces do not vanish when the final-state phase-space is suf- 
ficiently constrained. To deal with these terms we follow the prescription of Ref. [21]. We 
define the non-singlet axial current by removing 75 : 

tyfltfb -> — gj ^IvloTi^- (7) 

These expressions are equivalent in the <i — > 4 limit. Since the Levi-Civita tensor is a 
four-dimensional object, it must be combined with the matrix elements only after they are 
rendered finite. Although this seems to imply that computations should be performed with 
open Lorentz indices, which would be cumbersome in realistic calculations, this can be 
avoided. Since we are interested in the products of two traces which each contains a single 
75, we only obtain products of two Levi-Civita tensors of the form € ftl ^ 2fl3fl4 e VlU2l/31 ' 4 . These 
can be simplified using the identity 

e ^2^ eV1U2U3Vi = -det|^|, n = fii, ..(14, v = 1/1...Z/4, (8) 

which can be easily continued to d ^ 4. We write the determinant in an explicit form 
and contract the free indices with the matrix elements. The resulting expressions become 
functions of scalar products of the particle momenta. The replacement in Eq. © violates 
the Ward identity that relates the renormalization of the axial and vector currents, since this 
definition of 75 does not anti-commute with all of the 7^ 112 ill . The following additional 
finite renormalization of the axial contribution must be performed: 

Zr = 1-^C F + ^ (22C 2 - ^C F C A + 2 -C F N^ , (9) 



log(A) 
A 



(6) 
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where Nf denotes the number of active fermion flavors and Cf,a are the standard QCD 
Casimir invariants. 



The differential cross section for the partonic process qq — > Z + X — > e + e + X can be 
written as 

IT T ^ u 

\M\ 2 = (10) 

where q 2 is the invariant mass of the di-lepton pair and H^ u and L^ u are the hadronic and 
leptonic tensors, respectively. Since we can define infrared safe observables in QCD us- 
ing leptonic momenta, the lepton tensor L^ u is irrelevant for the cancellation of soft and 
collinear singularities. It is therefore natural to take the leptonic phase- space in four dimen- 
sions, rather than d dimensions, to simplify the calculation of matrix elements. However, we 
must take care when writing the scalar products of leptonic momenta with partonic momenta 
in the double-real emission corrections. We consider here the qq — ► ggZ — > gge + e~ process 
for illustration. The phase-space for the gge + e~ final state factorizes into the phase-space 
for qq — > ggZ and the phase-space for Z — > e + e~. We perform a Sudakov decomposition 
of the final-state momenta in terms of the incoming partonic momenta pi,p2- Denoting the 
momentum of one of the gluons as p%, we have 

P3 = a 3Pl + hP2 + PZT- (11) 

Because of momentum conservation, we must define one relative angle between p 3 T and 
Pzt to parameterize the ggZ phase-space and one relative angle between p e x and pzr to 
parameterize the Z — > e + e~ phase-space. In d = 4, each transverse phase-space is 2- 
dimensional, and these two angles determine the relative angle between p 3 T and p e T- In 
d ^ 4, the ggZ transverse phase-space is (d — 2)-dimensional; an additional angle <p is 
needed to define the relative orientation between the planes defined by Pzt,P3T and by 
Pzt, PeT- The scalar product must therefore be written as 

Pzt ■ PeT = PzTPeT (cos (j) 3Z cos (j) eZ + sin (p eZ sin 3Z cos 0) . (12) 

In the limit d — > 4, we have cos0 — > ±1, and this expression reduces to 
PzTPeT cos ((fi 3Z ± 4>ez)- This makes it explicit that given the two angles (f) 3Z and 4> eZ , the 
orientation of p 3T and p eT is completely determined in d = 4. 



After all three components of the hard-scattering cross section are combined, additional coun- 
terterms are needed to remove initial state collinear singularities. It is straightforward to extend 
the numerical approach described in [ 17] to obtain the desired results. 

We have essentially two checks on our calculation. First, considering different cuts on the 
leptonic transverse momenta and rapidities and on the missing energy for W production, we verify 
cancellation of the divergences in the production cross sections. Because the divergences start at 
1/e 4 at NNLO, the cancellation of all divergences through 1/e provides a stringent check on the 
calculation. We also check that the vector and axial contributions are separately finite, as required. 
A second check is obtained by integrating fully over the final state phase- space and comparing 
against known results for the inclusive cross section. We find excellent agreement with the results 
of |9|] for all partonic channels. 
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III. PHENOMENOLOGICAL RESULTS 



We have implemented our calculation into a numerical program, and we now discuss several 
phenomenological results obtained using this code. We first present the input parameters. We 
use the MRST parton distribution functions [22] at the appropriate order in a s . We use M z = 
91.1875 GeV, T z = 2.4952 GeV, Br(Z -> e+e") = 0.0336, M w = 80.451 GeV, T w = 
2.118 GeV,Br{W -> ev) = 0.1068. We set \V ud \ = 0.974, \V US \ = \V cd \ = 0.219, and |V CS | = 
0.996, and obtain \ V u b \ and \V c b\ from unitarity of the CKM matrix. We neglect contributions from 
the top quark; these have been shown to be small in the inclusive cross section For electroweak 
input parameters, we use sin 2 6w = 0.2216 and «qed(^z) = 1/128. We set the factorization and 
renormalization scales to a common value, \x r = /// = /i, and employ various choices of /x in our 
numerical study. To perform the numerical integration we use the CUBA package [ 23 ] . 



LHC 7+Z 



LHC 7+Z 
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FIG. 1 : The production cross section (left panel) and the acceptance (right panel) as functions of the lepton 
pseudorapidity cut for neutral current l + l~ production at the LHC. The two charged leptons are required to 
have pt > 25 GeV, and their invariant mass is constrained to 66 < M\+i- < 116 GeV. The dotted green 
lines refer to the LO result for fx = Mz/2 and \i = 2Mz, the solid red lines indicate the NLO result, and 
the dashed blue lines denote the NNLO result. We note that the \i = M z /2 and \i = 2M Z NNLO lines 
almost completely overlap and are nearly indistinguishable in both panels, and that both NNLO lines are 
completely contained within the NLO results. 



The identification of electroweak gauge bosons at hadron colliders typically requires cuts on 
the transverse momenta and pseudorapidities of the charged leptons, as well as on the missing 
energy for W^-boson production. For Z production, the invariant mass of the di-lepton pair is also 
restricted to suppress the importance of photon exchange. We first study the importance of the 
NNLO QCD effects for the cross sections and acceptances as a function of kinematic cuts for Z 
production at the LHC. In Fig.[H the neutral current l + l~ rate and acceptance at the LHC is studied 
as a function of a cut on the leptonic pseudorapidities. The NNLO results are absolutely stable 
with respect to scale variations, with residual uncertainties much less than one percent, and are 
completely contained within the NLO uncertainty bands. 

In Fig. El we present the neutral current l + l~ rate and acceptance at the LHC as a function of a 
minimum lepton px cut which we refer to as pr, c - Several comments regarding these results are 
required. 

• There is a kinematic boundary at p\ = M z / 2 above which the pure Z contribution to the LO 
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cross section vanishes in the limit Y z — > 0. At higher orders, soft gluon effects are important 
near this boundary. We expect the fixed-order result to be very accurate for values of pt, c 
away from this boundary. Evidence for this is provided by the close agreement between 



NLO and MC@NLO for a similar boundary at p b T = M w /2 in W production II ill . 



Below px,c — 40 GeV, the NNLO results are absolutely stable with respect to scale varia- 
tions, with residual uncertainties less than one percent, and are almost completely contained 
within the NLO uncertainty bands. 

For higher values of p Tc , there are large shifts when going from NLO to NNLO, and the 
scale uncertainties underestimate the corrections. This is not too surprising; in the limit 
Y z — > the LO result vanishes in this region since an additional radiated gluon is needed to 
have pt > M z / 2, and what we call NLO is the first term in the perturbative expansion. The 
absolute magnitude of the shift is also consistent with an O{o 2 s ) effect. 




Pt,= [GeV] p TiC [GeV] 



FIG. 2: The production cross section (left panel) and the acceptance (right panel) as functions of the lepton 
transverse momentum cut for neutral current l + l~ production at the LHC. The two charged leptons are 
required to have \rj\ < 2.5, and their invariant mass is constrained to 66 < Afj+j- < 116 GeV. The dotted 
green lines refer to the LO result for \i = Mz/2 and fi = 2Mz, the solid red lines indicate the NLO result, 
and the dashed blue lines denote the NNLO result. 

We now discuss a Tevatron analysis of the W -boson cross section. CDF recently presented a 
measurement of the W — > ev cross section in the forward rapidity region, 1.2 < \r\\ < 2.8, and 
compared this result to the central cross section [24]. Different values of Bjorken-a; contribute to 
each rapidity region, and measuring the central/forward cross section ratio may provide a useful 
constraint on parton distribution functions. The geometric and kinematic cuts in each region on 
the charged lepton pseudorapidity and transverse momentum, and on the missing energy, are listed 
below. 

• Forward: 1.2 < \rj\ < 2.8, E T > 20 GeV, E T > 25 GeV; 

• Central: \r}\ < 1.1, E T > 25 GeV, E T > 25 GeV. 

In the central cross section analysis there are additional selection cuts requiring the electron to be 
in the fiducial region of the calorimeter, and for the tracker to find an electron with p T > 10 GeV 
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consistent with the energy deposition in the calorimeter [ 12]. These cuts give an additional factor 
Acor = 0.6985, so that the acceptance in the central region is A cen = A geom x A cor . We compute 
A geom through NNLO in perturbative QCD and use the given A cor to determine the acceptance in 
the central region. 

We present in Table[Qthe predictions for the acceptances in the central and forward regions, and 
for the central/forward ratio R c /f. The uncertainties in this table have been obtained by computing 
the results for the scale choices M w /2 < fi < 2M W and equating the spread with the residual 
uncertainty. This procedure is supported by the NNLO results lying within the ranges indicated by 
the NLO scale variation. The NNLO theoretical uncertainties are at the 0.25% level or less, and 
are completely negligible. The forward region acceptance in particular is absolutely stable against 
radiative corrections. In the central region we observe an error reduction of a factor of five when 
the NNLO QCD effects are included. Only experimental errors and parton distribution function 
uncertainties remain, indicating that this measurement can potentially provide useful constraints 
on parton distribution functions. Our result Rfyj, = 0.9266(19) is in good agreement with the 
preliminary value obtained by CDF, R e * p f = 0.925(33) lET 





NLO 


NNLO 


A for 


0.2616(2) 


0.2614(2) 


A 


0.2458(28) 


0.2422(5) 


Rc/f 


0.940(12) 


0.9266(19) 



TABLE I: The theoretical predictions for the forward region acceptance Af or , the central region acceptance 
A cen , and the ratio of central/forward acceptances R c /f, together with their associated uncertainties at NLO 
and NNLO. We note that both the geometric acceptance and the factor A cor have been included in the 
central region result. 



IV. CONCLUSIONS 

In this paper we described a computation of the fully differential cross section through NNLO in 
QCD for W and Z boson production in hadronic collisions. Our result includes spin correlations, 
finite width effects, 7 — Z interference and allows for the application of arbitrary cuts on the 
final-state decay products. We have incorporated our result into a numerical code FEWZ available 
at the web site http://www.physics.hawaii.edu/~kirill/FEHiP.htm We believe this program will be 
invaluable for precision electroweak studies at both the Tevatron and the LHC. 

We studied several LHC and Tevatron examples where precise predictions for gauge boson 
acceptances are required. The theoretical prediction for neutral current / + / production at the 
LHC is absolutely stable with respect to residual scale dependence, with a remaining theoretical 
uncertainty much less than 1%, as long as the minimum pr cut on the leptons is less than the 
kinematic boundary value p b T = Mz/2. For momenta above p h T the LO result vanishes in the 
limit T z — > 0. The NNLO calculation provides the first radiative correction in this region, and a 
significant scale variation remains. 

We also studied the W cross section in the central and forward regions, recently analyzed by 
the CDF collaboration. The theoretical predictions for the acceptances in each region have residual 
uncertainties less than 0.25%. Our calculation of the ratio of central and forward cross sections, 
Rf = 0.9266(19), is in good agreement with the preliminary CDF result = 0.925(33). 
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